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Abstract: We calculate analytically the improvement coefficients of the static axial and 
vector currents in 0(a) improved lattice QCD at one- loop order of perturbation theory. The 
static quark is described by the hypercubic action, previously introduced in the literature 
in order to improve the signal-to-noise ratio of static observables. Within a Schrodinger 
Functional setup, we derive the Feynman rules of the hypercubic link in time-momentum 
representation. The improvement coefficients are obtained from on-shell correlators of the 
static axial and vector currents. As a by-product, we localise the minimum of the static 
self-energy as a function of the smearing parameters of the action at one-loop order and 
show that the perturbative minimum is close to its non-perturbative counterpart. 
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1. Introduction 



The hypercubic (HYP) link has been introduced a few years ago in order to study the 
static quark/anti-quark potential with improved statistical precision [1]. In this context it 
proved to be remarkably effective at large quark distances, where a description of the static 
propagators in terms of thin gauge links was known to lead to an excessive noise. The origin 
of the statistical improvement has been subsequently identified as due to a strong reduction 
of the static self-energy for appropriate choices of the HYP smearing parameters [2]. 

Since the binding energy of static mesons is amplified by self-energy contributions, the 
adoption of the HYP link as a parallel transporter in the static action, originally introduced 
in ref. [3], has allowed for lattice simulations with smaller binding potentials, thus triggering 
an exponential improvement of the signal-to-noise ratio of meson correlators at large time 
distances. This helped significantly in several applications of Heavy Quark Effective Theory 
(HQET) [4-9]. 

On a theoretical side, the HYP smearing procedure looks superior to other techniques 
in that it mixes gauge links within hypercubes attached to the original link only, which 
allows to preserve locality to a high degree. 

In spite of the statistical gain, the HYP link brings on an increase of operational com- 
plexity, which in some applications has compelled to resort to approximations. For instance, 
this is the case with the 0(a) improvement of the static-light bilinear operators, which has 
been addressed to in [2, 10]. In these papers the improvement of the static axial and vector 
currents is obtained by imposing that on-shell renormalised correlators of the currents with 
different choices of the HYP smearing parameters differ by 0(a 2 ) effects. Obviously, this 
procedure allows to relate the operator improvement coefficients corresponding to different 
static discretisations, but in order to isolate all of them, the improvement coefficients must 
be known for at least one choice of the smearing parameters. 

An approximate solution has been obtained by using a one-loop perturbative evalu- 
ation of the improvement coefficients with static discretisations characterised by simpler 
lattice Feynman rules, such as the original Eichten-Hill (EH) one or the so-called APE 
discretisation (see [2] and sect. 2 for definitions). The approximate solutions thus obtained 
for other choices of the HYP parameters are therefore affected by 0(<7q) systematic errors. 
As such, they have to be regarded as effective estimates, which are neither purely one-loop 
perturbative, nor fully non-perturbative. Of course, the amount of the 0(<7q) terms can be 
assessed only if the one-loop value of such improvement coefficients is known exactly. 

The aim of the present paper is precisely the analytical determination of the improve- 
ment coefficients of the static axial and vector currents at one-loop order of perturbation 
theory (PT) for the most convenient choices of the HYP link. In order to accomplish this, 
we adopt a formalism based on the Schrodinger Functional (SF) [12,13], derive the lattice 
Feynman rules of the HYP link at one-loop order of PT in time-momentum representation, 
and follow [10, 11] as for the choice of the operator improvement conditions. 

As a by-product, we also calculate the static self-energy at one-loop order as a function 
of the smearing parameters and perform a numerical minimisation of it. Remarkably, the 
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minimum of the self-energy is found for a choice of the smearing parameters which is very 
close to what is nowadays called the HYP2 action [2], thus showing that the self-energy is 
dominated by perturbation theory 1 . 

The paper is organised as follows. In sect. 2 we review some basic definitions and 
introduce our notation for the HYP static action. Sect. 3 is devoted to a discussion of the 
effects of the SU(3) projection, used in the HYP procedure, at one-loop order of PT; we 
show that some results, previously known in the literature, are largely independent of how 
the projection is defined. In sect. 4 we provide a simple derivation of the Feynman rules 
of the temporal HYP link in time-momentum representation. A systematic study of the 
one-loop static self-energy is described in sect. 5. In sect. 6 we collect our results for the 
improvement coefficients. Conclusions are drawn in sect. 7. For the sake of readability we 
strived to avoid as much as possible unessential technicalities. For this reason, algebraic 
details have been almost completely confined to the appendices. 

2. Notation and basic definitions 
2.1 HYP static action 

In order to set-up the notation, we review some basic formulae concerning the static theory 
and the HYP link. Static quarks are represented by a pair of fermion fields 
propagating forward and backward in time, respectively; their dynamics is governed by the 
lattice action [3] 

S stat [^ h ,^] = a 4 ^^ h (x)V^ h (x)-^(x)V ^(x)} , (2.1) 

X 

where the forward and backward covariant derivatives Vo, Vq are defined via 
VoV'h(^) = \ [W (x)^(x + a6) - ^(x)] , 

VoY> h (z) = \ [M*) ~ W °( x ~ aOrVhO* - o6)] . (2.2) 

The field V'h(V'h) can be thought of as the annihilator(creator) of a heavy quark. Simi- 
larly, V'h(V ; h) creates(annihilates) a heavy anti-quark. Each field is represented by a four- 
component Dirac vector, yet only half of the components play a dynamical role, owing to 
the static projection constraints 

P + V h = i> h , i^ h P+ = ? h , p+ = ^(i + 7o) ; 

P_^ E = Vh ! Y^-=^h, P=^(l-7o). (2.3) 

The static action has a functional dependence upon the parallel transporter Wo, which 
is assumed to be represented in this paper by a temporal HYP link. According to the 

lr The perturbative minimisation has been performed originally by R. Hoffmann using periodic boundary 
conditions, but it has never been published. We did it independently and are grateful to him for private 
communications. 
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original definition [1], the HYP link is obtained through a three-step recursive smearing 
procedure, where the original thin link is decorated with staples belonging to its surround- 
ing hypercube, i.e. 

W„(x) = W$\x) = P su(3) [(l - a l )U fl {x) + 

+ Y E WW(x)W$(x + d>)WW(z + aA)+] , (2.4) 

W^(x)=V su(3) [(l-a 2 )U^x)+ 

+ j E W$r(x)Wj&(x + ap)W$„(x + aii)i], (2.5) 

w£l p (x) = V svi3) l(l-a 3 )U fM (x) + 

+ Y Yl Ur,{x) Up(x + arj)U v {x + ap.?] , (2.6) 

where U-^{x) = uj jL {x — ap). In Eq. fl2,4|) denotes the fundamental gauge link and the 

(2) 

index v in W^-l indicates that the fat link at location x and direction p is not decorated with 
staples extending in direction v. The decorated link wf£l is then constructed in Eq. ( |2.5| ) 
with a modified APE blocking from another set of decorated links, where the indices pv 
indicate that the fat link W^ pv in direction p is not decorated with staples extending in 
the p or v directions. Finally, the decorated link wj^ pv is constructed in Eq. ( |2.6| ) from the 
original thin links with a modified APE blocking step where only the two staples orthogonal 
to p, v and p are used. After each smearing step, the new fat link is projected onto SU(3). 



Popular choices of the smearing parameters are represented by 

EH : ( ai ,a 2 ,a 3 ) = ( 0.0, 0.0, 0.0) , (2.7) 

HYP1 : (ai,a 2 ,Q3) = (0.75, 0.6, 0.3), (2.8) 

HYP 2 : (ai,a 2 ,a 3 ) = ( 1.0, 1.0, 0.5) . (2.9) 



The APE action of ref. [2] is characterised by a parallel transporter, which is not 
projected onto SU(3). Therefore, the action cannot be generated through a HYP link for 
any choice of the smearing parameters at a non-perturbative level. Nevertheless, the SU(3) 
projection is irrelevant at one-loop order of PT, as shown in next section. It follows that, 
as far as we are concerned here, the APE action is obtained from 

APE : (01,02,03) = (1.0, 0.0, 0.0) . (2.10) 

2.2 O(a) improvement of the static currents 

We adopt a SF topology where periodic boundary conditions (up to a phase 9 for the light 
quark fields) are set up on the spatial directions and Dirichlet boundary conditions are 



-3- 



imposed on time at xo = 0, T. We also assume no background field in our theory. The SF 
formalism has certain advantages, which have been widely discussed in the literature; we 
refer the reader to [14] for an introduction to the subject. In particular, it is here adopted 
since it allows for an easy on-shell formulation of the improvement problem. 

In this paper, we are interested in the static axial and vector currents 

Af at (x)=lp(x) 70l5 iJ h (x) , (2.11) 
y o stat (x)= ^ (x)7oV , h(x) . {2 .12) 

The 0(a) improvement of these operators with Wilson light and EH static quarks has been 
studied at one-loop order of PT in [10, 11]. The general improvement pattern is as follows: 

(A s R tat ) = Zf & \\ + 61 tat am q )(Af at ) , (2.13) 

(y R stat )o = ^ tat (l + 6^ tat am q )(Vi stat )o , (2.14) 

with 

(^ tat )o = < at + < tat o< at , 6A$*\x) = ?(*) 7 ,75^ (V, + V*) i> h (x) ; (2-15) 

(V^) = y stat + acT'SV^ , 5V stat (x) = ^(x) 7 ^ ( V, + V*) ^ h (x) . (2.16) 

The improvement coefficients c^ at , Cy at , 6^ at and 6^ at depend on the gauge coupling and 
are perturbatively expanded according to 

bT = \ + bT (l) gl + 0{gt) ; X = A, V . (2.18) 

We do not review explicitely the improvement conditions used in [10,11] here, but perform 
the same perturbative analysis with HYP action and refer the reader to those papers for 
details. In particular, we extract the one-loop improvement coefficients c^ at ^ and frj 3 ^ 1 ) Q f 
the static axial current as explained in appendix B of [11], and the one-loop improvement 
coefficients Cy at ^ and 6^ at ^ of the static vector current according to Eqs. (4.3)-(6.13) 
of [10]. 

3. SU(3) projection at one-loop order of PT 

The SU(3) projector "P S u(3) is n °t uniquely defined and the resulting HYP link depends 
in principle upon the chosen definition. Here we follow [2]: if A denotes a generic 3x3 
complex matrix, its projection M onto SU(3) is defined by 

M = ■ X , X = A(A^A)~ 1/2 , (3.1) 
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i.e. the matrix A is first unitarised, and then its determinant is rotated to one. Another very 
popular definition goes through the maximisation (minimisation) of the real (imaginary) 
part of ti(M^A). The latter has been employed in particular in [15], where two valuable 
propositions have been proved, which allow for a substantial simplification of perturbative 
calculations at one-loop order. The propositions concern the adjoint field of the HYP link, 
when this is expanded in power series of the adjoint (gluon) field of the fundamental gauge 
link. They read as follows: 

1. the linear term of the series is invariant under the action of Vsv(3)> 

2. the quadratic term is anti-symmetric in the gluon held. 

Although the proof proposed in [15] is based on the specific definition of the SU(3) 
projector as maxAf {Retr[MTyl]}, the above statements hold as well if Eq. (3.1) is adopted 
instead. A proof of this is straightforward. To fix the notation, we define the gauge field 
(x) of the fundamental link via 

U^x) = exp{a 50 <(z)T a } , (3.2) 

where {T a } a= i...8 denotes an anti-hermitian representation of the Gell-Mann matrices. 
Following [15], the unprojected APE-blocked links of Eqs. (|2.4|)-( |2"lf ) are written in the 
general form 

Vn{x) = 1 + ag Q ^f ljLU {x,y)q u (y) + a 2 gl ^ h^ p ( (y)q p (z) + 0(g$) , (3.3) 

u;y vp\yz 

with real vertices f,h depending on the level of the HYP smearing and q p (x) = q^(x)T a . 
It should be observed that, since V p is neither a SU(3) nor a U(3) matrix, no gauge field 
can be associated with it. Projecting the link according to Eq. (|3.1| ) leads to 

V„(x) \Vl{x)V p {x)\ = 1 + ago2^Uv(x,y)q„{y) 

a 2 2 



vp\yz 



,2„2 



and 



up;yz 

+ O( 5o 3 ) , (3.4) 

detjv^z) [Vj(^)^(x)]" 1/2 | = 1 + O( 5o 3 ) . (3.5) 

Eq. (|3.5| ) can been easily checked through an algebraic program, such as FORM [16]. The 
projected link Vsv(3)\Vp{x)\ is now in SU(3) by construction. Consequently, we are allowed 
to introduce a smeared gauge field 

?WW] = e W {agov a Jx)T a } . (3.6) 
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As a remark we also note that, since the perturbative expansion of the staples generates 
gauge fields at all orders in PT, v p cannot be assumed as a linear function of q^, and has 
to be perturbatively expanded in its turn, i.e. 



v + ag v" + 0{go) ■ 



Accordingly, Eq. (3.6) reads 



•,2„2 



Vsv (3) [V^x)} = 1 + a 90 < + a 2 glv'l + ^« + 0(gg) 



(3.7) 



(3.8) 



Equating Eq. ( |3.4| ) and Eq. (3.8) order by order in the bare coupling allows to express the 
smeared gauge field in terms of the unsmeared one. In particular 



O(g° ) : 1=1, 

O(ffo) : v 'l =\^Zh^ p {x,y,z)[q v {y),q p {z)} 

up;yz 



(3.9) 
(3.10) 
(3.11) 



Eqs. ( |3.10| )-( 3,11 ) correspond precisely to propositions 1 and 2 of [15]. It is worth noting 
that tadpole diagrams associated with Eq. ( 3. 1 1| ) do not contribute to any one-loop pertur- 
bative calculation in absence of a background field, although they are 0(<7q). Indeed, the 
commutator [T a ,T b ] = —f abc T c is anti-symmetric in a «-> b, while the Wick contraction 
(q a q b ) is symmetric. Nevertheless, one should avoid the conclusion that tadpoles do not 
contribute at all, because terms such as v'^ t (x)v' fl (x) have to be always considered. More- 
over, the vanishing of the above-mentioned contributions is not anymore true in presence 
of a background field, where PT becomes more cumbersome. 

To summarise, Eqs. (3.1C)-( 3.11| ) show that the smeared gauge field can be taken 
as a linear function of the unsmeared one at one-loop in PT and zero background field. 
The effect of the SU(3) projection can be disregarded in the linear term, and non-linear 
contributions should be discarded when deriving the Feynman rules of the HYP link. 



4. Feynman rules in time-momentum representation 

Translational invariance along time is broken within the SF formalism. This fact compli- 
cates perturbation theory, which is usually worked out in momentum space, and makes it 
natural to adopt a mixed approach, known as time-momentum representation, where only 
the spatial coordinates are Fourier transformed. Accordingly, tree-level propagators and 
interaction vertices are functions of time and spatial momenta. 

We derive the Feynman rules in time-momentum representation for the temporal com- 
ponent of the HYP link according to two different and independent procedures. The first 
derivation makes use of the results obtained in full momentum space [17] and arrives at 
time-momentum representation by inverse Fourier transform in time. This is possible in 
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presence of homogeneous Dirichlet boundary conditions, which allow to continue the SF 
periodically to all times with no non-trivial terms at the boundaries. The second derivation 
follows by direct construction upon representing the HYP smearing as a lattice differential 
operator acting linearly on the fundamental gauge fields. In this section we report only 
the first derivation; the second one, which is quite lengthy, is sketched in appendix A. For 
the sake of completeness, we report in appendix B the Feynman rules in time-momentum 
representation for the spatial components of the HYP link, which may be useful for some 
applications of HQET (e.g. see ref. [5]) and for dynamical n-HYP [18]. 

In order to fix the notation, we define the smeared gauge fields of the three levels of 
the HYP procedure according to 



wW(x) = exp{a 50 4 fe) (x)}, 



k = 1,2,3 . 



(4.1) 



Under the assumption of periodic boundary conditions in space and homogenous Dirich- 
let ones in time without nontrivial boundary terms in the action, the gauge field of the 
fundamental (and smeared) link admits a 4-dimensional Fourier transform 



v 



— < Pfj, 
a 



TT 



2ir 
L a 



(4.2) 



where L is the toroidal extension of the space-time dimensions and the additional phase 
shift in direction fi is due to the very definition of the gauge fields, which are supposed to live 
between neighboring lattice points. Factorising the loop-sums over the spatial directions 
allows to express the gauge field in time-momentum representation in terms of the fully 
Fourier-transformed one, i.e. 



PO 



F3 5^e ipx g (zo;p) , 



L 3 



Qk{x) = ^3 E e ipx e*f j £ eiP ° X °^(p) = ^3 E e ipx e^g fc (x ; p) . 



L 3 

P L Po 

By direct comparison, it follows 



9o(xo;p) = ^E eW ° elfpo 9oW . 
%(x ;p) = ^E eW0 ^W • 



(4.3) 
(4.4) 



(4.5) 
(4.6) 



We now consider the Feynman rules in momentum representation [17], 

¥(p) = &f(*(p) + O(flo) , 



(4.7) 
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where 



(4- 



^(p) = l + a 2 (l + a 3 )-f (l + 2a 3 )a 2 (p 2 -pl-pl) + ^ ]J a 2 p 2 . (4.9) 



i (3) 

Since we are interested in the static propagator, we focus on 5g , for which we get 



po 



PO 



foo(p)qo(p) + ^2 fok(p)Qk(p) 



k=l 



(4.10) 



We then observe that 

3 

/oo(p) = 1 - y ^2( a2 Pk)^ok(p) , 



(4.11) 



k=l 



Ct2 



n ok (p) = 1 + a 2 (l + a 3 ) - y(l + 2a 3 )(a 2 p 2 + a 2 ]5 2 ) + 



, a 2«3, 2~2\/ 2~2\ 

+ —^( a PiA a Pi) 



kkip) = ^r(aPo)(aPk)tt ok (p) 



0,i) ik 



(4.12) 
(4.13) 



The above equations show in particular that figfe an d /oo depend only upon the spatial 
components of the Fourier momentum. Therefore, 

^S 3) (^o;p) = /oo(p)<?oOro;p) + 

3 

+ y E( a ft) fl ^P)7E eWel> (aPo)Up) + O( 50 ) • (4.14) 



k=i 



Po 



The second term on the right hand side of the previous equation can be easily worked out, 
and we finally obtain 



Bq\x ; p) = /oo(p)<?o(zo; p) - ^ ^2(ap k )n ok (p)d qk{xo; p) . (4.15) 



iot\ 
~6 



fe=i 



A very compact way of writing the above expression is by introducing an effective vertex 
and auxiliary indices, i.e. 



-Bq 3) (^o;p) = ^2Vo- i (p)q tl ( i) (x + as{i);p) , 



(4.16) 



i=0 



all collected in Table |l|. Comparing the latter with Table 4 of [2] shows that the Feynman 
rules of the HYP link resemble closely those of the APE one, from which they differ by the 
presence of a form factor f^ofc- 
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i nii) s(i) Vb;i(p) 

o o o i-tEli« 2 p 2 fc ^(p) 

1,2,3 i + if a p ti{i) n 0fl{t) (p) 

4,5,6 i-3 1 - 1 f 1 op jU ( i )O 0A1 ( i )(p) 

Table 1: Fcynman rules of the temporal HYP link in time-momentum representation. 

5. The static self-energy at one-loop order in perturbation theory 

The binding energy E stSu t of a static-light meson controls the exponential decay rate of 
the associated two-point function. It depends upon the choice of the parallel transporter 
and diverges linearly in the continuum limit. As such, it can be perturbatively expanded 
according to 

£stat ~ E 8e u + 0(a°) ~ -e^g 2 + . . . . (5.1) 

a 

The coefficient represents an ultraviolet property of the static action. Accordingly, it 
is insensitive to the specific correlation function from which it is measured, as well as to 
the choice of the boundary conditions. 

In appendix A of ref. [2], e^ 1 ' is defined from the SF boundary-to-boundary correlator 
f^ h , made of two static quarks propagating across the bulk region. Here we adopt a 
different definition and extract from the boundary-to-boundary correlator /f tat , made 
of one static and one relativistic propagator 2 , i.e. 

/•stat(l) 

e (D = _ Hm . (5.2) 

where fi*^ denotes the sum of the two Feynman diagrams depicted in Fig. [I], with the 
light-quark line kept at tree-level. Evaluating Eq. Q5.2j ) leads to the expression 

^1^2^3=0 

with non-zero coefficients e^^ fc collected in Table ^. A derivation of Eq. ( |5.3| ) is reported 
in appendix]^. The static self-energy is a multivariate polynomial of (ai, 02,0(3). The 
structure of the polynomial is triangular, i.e. the only non-vanishing contributions are at 
< &3 < k2 < k\ < 2. This is an obvious consequence of how the three steps of the HYP 
smearing are defined in Eqs. Q2.4| )— ([2-6D . 

It is also worth noting that higher order perturbative corrections to Eq. Q5.1| ) have the 
same functional form (but different coefficients) as Eq. (|5.3]), since no HYP links appear 



2 For the original definition of /f a \ see Eq. (3.23) of [11]. 
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Figure 1: Feynman diagrams contributing to the static self-energy at one-loop order of PT. Eu- 
clidean time goes from left to right. Single (double) lines represent relativistic (static) valence 
quarks. 



in the light quark action, and therefore dynamical quark loops are directly connected to 
thin gluons only. This is not true anymore if the HYP link (or variants of it) enters the 
covariant derivatives of the Wilson action, such as in ref. [18]. In this case, dynamical 
quark loops may have a functional impact on the static self-energy. 

A numerical minimisation of the coefficient as a function of the smearing parame- 
ters is easily performed. To this aim we used a MATLAB routine. The global minimum is 
found for a choice (a*, a^, a^) given by 

<£L = e (1) | K , a ., Q *) = 0.03520(1) , (aj.o&og) = (1.0, 0.9011, 0.5196) . (5.4) 

A two-dimensional contour plot of e^ 1 ) as a function of ai and «3 at a.\ = 1.0 is 
shown in Fig. |2[ Remarkably, the perturbative minimum is obtained with a choice of the 
smearing parameters which is very close to the HYP2 definition, see Eq. ( |2.9P , found via a 
numerical minimisation performed at a non-perturbative level [2,19]. The closeness of the 
perturbative and non-perturbative minima suggests that the static self-energy is dominated 
by PT, as intuitively expected. A perturbative comparison between Eq. ( |5.4[) and HYP2 is 
also possible: we find Chyp2 = 0.03544(1), which is slightly higher than the result reported 
in Table 1 of [2] and differs by 0.7% from the perturbative minumum. 

6. Determination of the improvement coefficients 

Once the Feynman rules of the HYP link are known in time-momentum representation, the 
extraction of the improvement coefficients of the static axial and vector currents at one-loop 
order of PT can be performed along the lines of refs. [10, 11]. As already warned in sect. 2, 
we do not review here the improvement conditions (and related perturbative equations) 
used in those papers, since our analysis has nothing new to add from a methodological point 
of view 3 . Instead, the reader is referred to appendix B of [11] as regards the improvement 
coefficients c^ at ^ and b^^ 1 ^, and to Eqs. (4.3)-(6.13) of [10] as for Cy at ^ 1 '' and by at . 

However, we observe that, while reproducing the calculation of b A with EH action, we have found 
a little mistake in Eq. (B.20) of [11]. This equation should be replaced by 

= (1 + &r (0) am( 0) ) [x W + Z^ t (1) X (0) ] + bT t(0) am^X m . (6.1) 

and the correct value of the improvement coefficient is 5j at W _ 0.0041(4) with EH action and 5j at W _ 
0.0925(6) with APE action. Note that the latter differs from Eq. (3.8) of [2], which has been determined 
using the wrong value of &l tat(1) [EH] as an input. Note also that all these values for j^ 3 ^ 1 ) are numerically 
rather small. The mistake has been discussed with and recognised by the authors of refs. [2, 11]. 
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a3 0.5 
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a2 

Figure 2: Countour plot of the coefficient as a function of a2 and a?, at ai = 1.0. The colour 
spectrum goes from violet (higher values) to red (lower values). 

Upon this premises, we collect our results in Table |||. It is understood that all the 
conditions employed define the improvement coefficients up to O(a) terms and need to be 
extrapolated to the continuum limit. 

For what concerns the axial current, we obtain the coefficient from three inde- 

pendent improvement conditions, i.e. at € {0.5,0.75,1.0} and renormalised quark mass 

statfl) statfl) 

z = Lm R = 0.0. Having determined c A , we obtain b A from three new indepen- 
dent conditions, i.e at 9 G {0.0, 0.5, 1.0} 4 and z = 0.24. Continuum approach is shown in 
Figs. |-| 

As for the vector current, we obtain the coefficient cv at ^ from three different im- 
plementations of the axial Ward Identity, corresponding to pairs of 0-angles (61,02) £ 
{(0.0, 0.5), (0.0, 1.0), (0.5, 1.0)}, and the coefficient fr^ at W f rom the ratio of a three-point 
correlator of the static vector current at renormalised quark masses (#1,22) = (0.0,0.24) 
and three values of the #-angle, viz. 9 G {0.0,0.5,1.0}. Convergence to the continuum 
limit is shown in Figs. ||-|| 

In all cases, different definitions lead to a consistent continuum limit. The uncer- 
tainty on the final numbers has been estimated as the maximal difference of the continuum 
extrapolations corresponding to independent definitions. Data fits have been performed 
according to the asymptotic expansion 

c (!) - c + * (!) + * (!) * (!) + * (i) 2 + B > (!) 2 * (!) + • ■ ■ <"> 

with c = c f t(1) , &f t(1) , 4 tat(1) , C t(1) ■ 



4 It is not possible to obtain at 9 = 0, since the operator <5^4o tat vanishes at tree-level for this 

particular choice of the (9-angle. 
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[hJ, k \ 

[0,0,0] 
[1,0,0] 
[1,1,0] 

[1,1,1] 
[2,0,0] 
[2,2,0] 
[2,2,1] 
[2,2,2] 



e (1) 

0.168487(1 
-0.222222(1 
-0.041164(1 
-0.015484(1 

0.111111(1 

0.023521(1 
-0.002620(1 

0.019055(1 



action 


v 
A 


stat(l) 


,stat(l) 
Ox 


HYP1 


A 


0.0029(2) 


0.0906(2) 


HYP1 


V 


0.0223(6) 


-0.0212(8) 


HYP 2 


A 


0.0518(2) 


0.142(1) 


HYP 2 


V 


0.0380(6) 


-0.0462(8) 



Table 2: Coefficients of the static self- 
energy at one-loop order of PT. 



Table 3: Improvement coefficients of the 
static-light currents at one-loop order of PT. 



Since O(a) lattice artefacts depend upon all the details of the calculation, the infor- 
mation given in the previous paragraphs has to be complemented with further technical 
details, in order to allow for complete reproducibility. In particular, plots reported in 
Figs. refer to a choice of the critical quark mass as obtained from the PCAC relation. 
Numerical values used have been taken from refs. [20,21]. Moreover, the determination 
of 5^ tat ( 1 ) requires the one-loop coefficient z^ at Q f the renormalisation constant of the 
static axial current as input. The scheme dependence of this coefficient has no effect on 
the continuum limit of , but it changes its continuum approach. Here, for each sin- 

gle 0- value, we used a definition of ^ from Eq. (4.14) of ref. [11], where the 0(a 2 ) 
^-dependent corrections at finite lattice spacing are taken as part of the definition. This 
allows for a strong cancellation of the O(a) lattice artefacts in ^ tat ( 1 ^ an d leads to a safe 
continuum extrapolation. Finally, the particular implementation of the axial Ward Identity 
needed for cv at ^ , which has been adopted here, refers to a choice of the lattice topology 
as T = 2, according to the notation of ref. [10]. 

Once the improvement coefficients have been determined, one can look at the residual 
0(a 2 ) lattice artefacts of the observables from which those coefficients are extracted. As 
an example, we consider the cutoff effects 

Y^ 1 ( n I T \ f stat _i_ „stat f stat 

5XUa/L) = ±^q±l _ i xl = Z! tat Ja + ^J^ , (6.3) 
RV 1 ' X^(0) y/jf^ 

of the ratio in the chiral limit. Here we choose to renormalise the static axial current 
in the "lat" scheme, i.e. with Z s £ &t having only divergent logarithms without finite parts. 
Eq. (|6.3| ) can be expandend in PT, 

i; = ir+W+0( 3o 4 ), (6.4) 

5X1 = SX^ + gl5X l K w + O( 5o 4 ) , (6.5) 
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Figure 3: Continuum extrapolation of the improvement coefficient c^ at ^ with HYPf (left) and 
HYP2 (right) static discretisation. On each plot, the three curves refer to independent determina- 
tions at 9 — 0.5 (circles), 9 = 0.75 (triangles) and 9 = 1.0 (squares). 
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Figure 4: Continuum extrapolation of the improvement coefficient ^ tat ^ 1 ) with HYP1 (left) and 
HYP2 (right) static discretisation. On each plot, the three curves refer to independent determina- 
tions at 9 = 0.0 (circles), 9 = 0.5 (triangles) and 9 = 1.0 (squares). 



6X^ m 



Xk m (a/L) 



X™ (0) 



(a/L) X l R w (0) 
^ (0) (a/L) XT (0) 



X-, 



1,(0) 



a/L) 



X^ (0) 



(6.6) 



(6.7) 



We note that 5X^ 1 ' does not depend upon the choice of the static regularisation. In 
order to compare lattice artefacts corresponding to different actions, one has to consider 
at least the one-loop contribution 8XT\ which is plotted in Fig. for EH and HYP2 at 
9 € {0.5, 0.75}. Remarkably, 0(a 2 ) lattice artefacts are less than 0.1% for both actions. Al- 
though in principle the smearing of the gauge link could be responsible for an enhancement 
of the cutoff effects, we do not observe any sign of this in our observables. 
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Figure 5: Continuum extrapolation of the improvement coefficient Cv with HYPf (left) and 
HYP2 (right) static discretisation. On each plot, the three curves refer to independent determi- 
nations at (6>i,6>2) = (0.0,0.5) (circles), (01,0a) = (0.0,1.0) (triangles) and (0i,0 2 ) = (0.5,1.0) 
(squares). 




Figure 6: Continuum extrapolation of the improvement coefficient 5^ tat ( 1 ) with HYP1 (left) and 
HYP2 (right) static discretisation. On each plot, the three curves refer to independent determina- 
tions at 9 — 0.0 (circles), 9 = 0.5 (triangles) and 9 — 1.0 (squares). 

It is also important to stress that the perturbative values collected in Table |3] need 
not to be in agreement with the findings of refs. [2, 10], where a hybrid technique has 
been adopted, which mixes one-loop perturbative inputs obtained from the EH or APE 
discretisations with non-perturbative simulations of HYP static fermions. As an example 
of this, we consider the case of 6y at , for which the mixed procedure gives 

Cp! ~ \ ~ 0-014(3) 5o 2 + O( 5o 4 ) , (6.8) 

Cp 2 ~\~ 0-096(8) 5o 2 + O(g 4 ) . (6.9) 
As it can be observed, the exact perturbative coefficients are sensibly different from the 



-14- 



- -0.0004 



th = 
th = C 


0.5, EH 
.5, HYP2 


1 1 1 1 


1 1 1 1 


i i i i . 


A 






V A 

1A A £ 

o 

c 


A 


A 








o 



















.... 


.... 


.... 






0.002 0.004 0.006 0.008 0.01 

(a/L) 2 



-g 0.0002 



th = 
th 


75 EH 


i i i i 


i i i i 




i 7c LiVDn 
L/Oj n ird 


o 








£ 




A 


A 
















.... 


.... 


.... 





0.002 0.004 0.006 0.008 0.01 

(a/L) 2 



Figure 7: Comparison of the residual 0(a 2 ) lattice artefacts of the improved ratio at one- 

loop order of PT in the chiral limit, with 9 = 0.5 (left) and 9 = 0.75 (right) for the EH (circles) 
and HYP2 (triangles) actions. 

effective ones, which signals the presence of non-negligible 0(<7q) terms within the latter, 
quantifiable as the differences among the perturbative and the effective values. The hybrid 
procedure is such that the differences absorb 0(<7g) terms of both actions used in the 
improvement condition, for the limited range of the bare coupling (6.0 < (3 < 6.5) where this 
has been implemented. A larger discrepancy between the exact improvement coefficients 
of the two discretisations at 0(<?q) induces a larger discrepancy between the perturbative 
coefficients of Table |3| and their effective partners. However, we remark that the impact of 
these 0(<7q) terms is suppressed when 6y at is multiplied by a reasonably small quark mass 
arriq. Similar considerations can be done also for the other improvement coefficients. 



7. Conclusions 

In this work, we have studied the O(a) improvement of the static axial and vector currents 
at one-loop order of perturbation theory with hypercubic static and Wilson light quarks. 
Our methodology for the on-shell O(a) improvement is based on the Schrodinger Func- 
tional. The calculation is useful in that it allows to quantify the amount of the 0(<?q) 
terms present in the hybrid perturbative/non-perturbative procedure of refs. [2,10]. 

The lack of knowledge of the Feynman rules for the hypercubic static propagator 
in time-momentum representation has prevented such calculation over the past years. We 
simply observed that, since the static propagator is a product of temporal hypercubic links, 
which are not affected by the boundaries of the Schrodinger Functional, the Feynman rules 
need not to be necessarily derived from scratch, but can be obtained by inverse Fourier 
transforming in time the ones obtained in full momentum space, the latter being known 
from the literature [17]. 

Aside the improvement coefficients, we have found an analytical expression for the 
static self-energy as a function of the smearing parameters at one-loop order of perturbation 
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theory. This expression has been minimised through an optimal choice of the smearing 
parameters. We have found that the perturbative minimum is very close to the non- 
perturbative one, obtained by means of an analogous numerical procedure performed at 
non-perturbative values of the bare coupling [2, 19]. This is interpreted as a hint of fast 
perturbative convergence. 
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A. Derivation of the temporal Feynman rules by direct construction 

Owing to the irrelevance of the SU(3) projection at one-loop order of PT, the HYP smearing 
procedure can be considered as a three-step linear lattice differential operator acting on 
the gauge fields at subsequent levels. The perturbative expansion of Eqs. (O)-( p^ ) reads 

B&>(x) = (1 - a l )q,(x) + ^Y1 + 

+ Y E (x) - dld,B^{x)} + O( 50 ) , (A.l) 

B$,(x) = (1 - a 2 )q,(x) + f E fl $p(*) + 

+ jE a\d*d p B%(x) - $0^gU*)} + O(go) , (A.2) 

B%{x) = q,{x) + y E a2 K^(*) - ^q v (x)} + O(g ) , (A.3) 

where d p and d* denote respectively the standard forward and backward lattice derivatives. 

Since in the SF perturbative calculations are naturally performed in time-momentum 
representation, the above formulae must be Fourier-transformed along the spatial direc- 
tions. Before any explicit calculation, it is worth first deciding what is relevant to our aims. 

~ (3) I 1 

To be concrete, we are interested in B\ . Accordingly, the Lorentz sums in Eq. ( |A.lj ) are 

(2) (2) (2) 

purely spatial. The only components of B^-l, to be considered are B^,i and B^ with 
k = 1, 2, 3. Analogously it can be argued of the following steps. If the notation 11 A ~~> B" 
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means that B is needed to calculate A, then the whole calculation is summarised by 

^S 3) ~» Sgj, A; = 1,2,3; (A.4) 

5g - <L = 1,2,3; (A.5) 

<o - Bgi, fc^/ = 1,2,3. (A.6) 

In order to obtain Bq we first consider the Fourier transforms of Eqs. (A.1)-( [A3| ) 
separately. Afterwards, we insert the expression obtained for the first-level smeared gauge 
fields Bj^J, p into B^J, and the latter into B^ . The Fourier transform of the first-level 
smeared gauge field is given by 

•^0;^( x o;p) = 1-f a2 P™ koOo;p) (ap m )ad <Zm(zo;p) , (A.7) 

^ m=£U ' m^kl 

^fc-cU z o;p) = \ 1 - y f*( x o;p) + Y ( a Pm)( a P*0<7m0co;p) , (A.8) 

^ rrvfkt ' m^kl 

^ofc^o; p) = \ i - y f^^o; p) + y ( a ^( a P/)<?m(x ; p) , (A.9) 

where = (2/a) sin(apfc/2) denotes a lattice momentum. In particular, it should be 
remarked that m ^ k,£, i.e. the sum over m contains just one term. Similarly, one finds 

Bgl(x ; p) = (1 - a 2 )g (x ; p) + y E^ 2 " ^^SlK P) + 



4 

i^k 



i— 2^(ape)adoBy h {x ; p) , (A.10) 



-^i 2 o( x o; p) = (i - a2)%(^o; p) + y X]( 2 ~ a2 P 2 )^i!d<?( x o; p) + 

e^k 

+ 7Lww B /!0fc( a; i)iP) ' ( A - n ) 

tj^k 

as for the second-level smeared gauge field. Finally, 

3 

B^ixo; p) = (1 - ai)9o(a:o; p) + yE( 2 " a2 Pk)^l(x ; p) + 

fe=l 



3 

i^J^(ap fc )a^,S^(a;o;p) . (A.12) 



6 



Instead of working out the whole expression alltogether, we observe that the final result 
is expected to be a linear combination of qo(xo;p), 9fc(^o;p) and qk(%o + &;p). Single 
contributions can be considered separately. 
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In order to extract the coefficient of qo(xo,p), we observe that B^L and B^). have no 
temporal gauge field, which appears only in B^ kt Accordingly, B^/ Q is independent of go- 

~ (2) (2) ~ (2) 

Instead, B^. k has such dependence. The coefficient c multiplying q~Q in B^. k can be easily 
isolated, i.e. 



o 



1--T- (l + a 3 )2^ a Pi + —7- 11 a P€ • ( A - 13 ) 



4 V " 4 



Analogously, the coefficient c multiplying go m ^o can De isolated upon replacing B . k 
with its explicit value. Some algebra leads to 

4 3) = (l - 01 ) + £ ^(2 - a a #)4 2 > = /oo(p) • (A.14) 



k 

Mi 



6 

A similar procedure allows to extract the spatial contributions to Bq J ' (xq ; p) . 
B. Spatial parallel transporter 

In this appendix we report the Feynman rules for the spatial HYP link in time-momentum 
representation. We follow sect. 4 and base our derivation on the inverse Fourier transform of 
the Feynman rules in momentum space, first obtained in ref. [17]. We start from Eq. (fO), 

B { k ] (p) = JMUP) + E fkv(p)UP) + fkk(p)q k (p) + 0(fl&) • (B.l) 

v¥=o,k 

Contributions on the right hand side are worked out separately. Some algebra leads to the 
expressions 



^ 3;1) (xo;p) = ~£e™*°/ko(p)gb(p) 

PO 



-i—(ap k )^lko{p)adQqo(xo;p) , (B.2) 



= (ap k ) £ [Ag(p) + Ag(p)a 2 9 *a ] p) , (B.3) 

Bf 3) (x ;p) = ^e^'A fc (p)%(p) = 



l + y^o(p)a 2 3 *9 - 



^ (ap„) [Ag(p) + Ag(p)a 2 a *9 ] p) , (B.4) 



v¥=o,k 
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i fi(i) s(i) 





1 

2 k 

3 k 

4 k 

5 rj) Ok 

6 7]; Ok 

7 77; OA; 







l-% 



^fifco(p) + 



^ap fc O fc o(p) 
|Lap fc O fc0 (p) 

W«^[2Ag( P )-Ag(p)] 



iHfco(p) - E^ ,fc °&j a 8(p. 







1 



-1 



4o(p) -E^o,fc a ^ A £(p) 
ap fe [Ai; ) (p)-2Ag(p)] 

opfcA^(p) 

aj5fcAg(p) 



Table 4: Feynman rules of the spatial HYP link in time-momentum representation. The symbol 
77; Ofc denotes a sum over r\ ^ 0, k. 



where we have introduced the symbols 



A S(P) = -f ( a Pv) 1 + "2(1 + a 3 ) - + 2a 3 )(a 2 ^ fcir) ) 



Ag(p) = Tr(aA,) -r(i + 2a 8 ) 



Eqs.(| 



can be written in the compact form 

6(3) 



■ 1 2 -2 



k,r]) 



(B.5) 
(B.6) 



-Bfc Oo!P) = ^2 V k;i{v)%(i){ x o +as{i);p) 



(B.7) 



i=0 



where the vertex Vfc ; j(p) and the auxiliary indices /j, s are collected in Table |j. 



C. Derivation of Eq. ( gj ) 

To conclude, we detail the calculation of the static self-energy and provide a derivation of 
Eq. ( |5.3| ). We consider the boundary-to-boundary correlator /f tat , defined in ref. [11] via 



,12 



ft 



stat 



2L 6 



(C / (u) 7 5CUv)Ch(y)75C( Z )> 



(C.l) 



uvyz 



Expanding J| tat at one-loop order of PT allows to write the coefficient see Eq. (|5.2j), 
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in terms of the Feynman diagrams of Fig. [1[ i.e. 

a T wo 6 



x ^(i)/x(j)(^o -a + as(i),v -a + as(j);p) 

T wo 



/i 111 ^73 X] X] & ( n o,«o)V(n ,wo,p) , (C.2) 

/ L — >0 o 1j 



where d llu {xQ,yo,p) denotes the gluon propagator in time-momentum representation, cf. 
ref. [13] for a definition. The weight-coefficient b(uo,vo) is given by 

f 1/2 u = v , 
6K^o) = < (C3) 
[ 1 otherwise , 

and the interaction blob 

6 

V(u ,v ,p)= ^ S^^Vo ; i(p)Vo ; j(-p)d^^)(ii^ -a + as(i),v -a + as(j);p) (C.4) 

i,j=0 

denotes a HYP gluon propagating on the lattice from time no to time vq with spatial 
momentum p. The above expression can be simplified by using spatial rotational invariance, 
i.e. d\\ = c?22 = ^33 and dij = if i ^ j. Accordingly, the vertex reads 

V(u ,v ,p) = |Vb;o(p)| 2 doo(>o - a, v - o;p) + 

+ [l^0;l(p)| 2 + !^0;2(P)| 2 + |Vb ;3 (p)| 2 ] X 

x [dkk(u - a, v - a;p) + d kk (u ,v ;p) - 

- d kk (u ,v - a;p) - d kk (u - a,v ;p)] ■ (C.5) 

From Eq. ( p.5| ), we conclude that the HYP vertex Voy'(p) enters the coefficient only 
in the rotationally symmetric combinations 

h {t) (p) = \V 0;0 (p)\ 2 , 

h (s) (p) = \V 0; i(p)\ 2 + |Vb ;2 (p)[ 2 + |Vb ;3 (p)| 2 • (C6) 
These quantities are multivariate polynomials of the HYP smearing parameters, i.e. 

~ kl ik=o ^ ' ^ ~ ** !A * ! ** ! dcfdafdafr ' ( } 

(t s) 

The coefficients w k ^ k2ks have been written according to the generic Taylor expansion in 
several variables. They can be algebraically evaluated; the non-vanishing ones are reported 
in sect. C.l in units of the lattice spacing. Spatial rotational invariance is evident. 
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Upon inserting Eq. (|07|) into Eq. (^J) and Eq. ( |05|) into Eq. (|(l2D , the final result 
represented by Eq. ( |5,3| ) is obtained with coefficients 



4 a 



T u 



x { wi Sk 2 k 3 (P) d oo(u ~a,v - a;p) +} 

+ ^^(p) t dfcfc ( u o ~ a ' w o - a; P) + dfcfc(«o, vo; p) 
-d kk (u ,v - a;p) - 4fc(«o - a,v ;p)]} . 



(C.8) 
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C.l Coefficients w k t f 2k3 

4oo = 1 , (C9) 

™^ = -\{pl + pl+vl), (cio) 
«iio = -\(p\+pI + pD + \ (p\pI + p\p\ + p\p\) , (en) 
wfA = -\{pI+pI + pI) + \ (pIpI + pIpI + pIpD - - A vlvlvl , (C.12) 

*4oo = ^ (Pi + Pt + P3) + ^ (p?pl + pfpl + pip§) , (C.13) 
w 2io = Y^iPi + P2 + Pi) + \ {PiPl + f\P\ + pIpD ~ 

- ^ (p!p1 + pIp! + p!p1 + pIpI + ^p? + pip 2 ,) - 1 P ^^2 ) (ai4) 

4?i = ^ (pi + pt + + 1 (m + ^ + m) - 

- (pip! + pIpI + P2P1 + P2P3 + P3P1 + P3P2) - g P1P2P3 + 

+ Ya ( ^'p3 + pIpipI + pIpIpI) , (C15) 

w 220 = j^j (Pi + P2 + P3) + ^(PiPs + P1P3 + P2P3) - j^PipIpI - 

- gg (P1P2 + P1P3 + P2P1 + P2P3 + P3P1 + P3P2) + 

+ Y 2 ^\vIpI + P2P?P3 + pIpIpI) + ^ (pfpl + pfpl + pIpI) , (C16) 

™221 = ^ (Pi + P2 + P3) + 9 (PlPl + PlP3 + P2P3) - 

- ^ (P1P2 + P1P3 + P2P1 + P2P3 + P3P1 + P3P2) - \ p\p\p\ + 
+ ^ (P1P2 + P1P3 + P2P3) + ^ (P1P2P3 + P2P1P3 + P3P1P2) ~ 
-^(pfpfp|+pfp|pI+P»?), (C17) 

W 222 = (Pi + P2 + P3) + jg (PiPl + PlP3 + P2P3) - 

- ^ p\p\p\ - (pIpI + P1P3 + P2P1 + P2P3 + P3P1 + P3P2) + 
+ ^ (P1P2 + P1P3 + P2P3) + ^2 (PiPzPs + P2P1P3 + P3P1P2) - 

- ya {p\pipl + p\Pipl+pipip\) + ^pIpIpI . (c.18) 



22 



(s) 
^200 


— (p\ +P9 + > 


(C.19) 


(s) 

Win = 
Z1U 


7: (Pi + P\ + Pi) ~ 777 (PiPa + P1P3 + pIpI) , 


(C.20) 


(s) 
W 211 — 


g IPl + P2 + P3J g \PlP2 + Pl^3 + P2P?,) + ^2 PlP2P3 ' 


(C 91 


(s) 
^220 ~~ 


i (p? + pi + vl) - i (p?p| + p?p| + p|p|) + 




+ 


(P1P2 + P1P3 + P2P1 + P2P3 + P3P1 + P3P2) + 7^ P1P2P3 . 


(C.22) 


(s) 
™221 = 


1 (Pi +Pl+Pl)-\ (PiPl + PiPl + pIpD + 





+ ^ (p\p\ + P\P\ + P2P1 + P\P\ + P3P1 + P3P2) + 3 P\P\p\ - 

- ^ (P\p\p\ + P2P1P3 + v\Pip\) , (C23) 

w 222 = \(Pi+Pl+Pl)-\ (PiPl + PiPl + PiPl) + 

+ (P1P2 + P1P3 + P2P1 + P\P% + P\Pi + P3P2) + 3 P\P\p\ + 

- \ (P\P\P\ + PiPlPl + PiPiPl) + 4 (P\P\P\ + P\v\P\ + • (C24) 
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